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ABSTRACT 

Simplified  equations  for  slender  jets  are  derived  for  a  circular  jet  of  one  fluid  flowing  into 
an  ambient  second  fluid,  the  flow  being  confined  in  a  circular  tank.  Invistid  flows  are  studied 
which  include  both  surface  tension  effects  and  Kelvin-Helmholtz  instability.  For  slender  jets 
a  coupled  nonlinear  system  of  equations  is  found  for  the  jet  shape  and  the  axial  velocity 
jump  across  it.  The  equations  can  break  down  after  a  finite  time  and  similarity  solutions 
are  constructed,  and  studied  analytically  and  numerically.  The  break-ups  found  pertain  to 
the  jet  pinching  after  a  finite  time,  without  violation  of  the  slender  jet  ansatz.  The  system 
is  conservative  and  admissible  singular  solutions  are  those  which  conserve  the  total  energy, 
mass  and  momentum.  Such  solutions  are  constructed  analytically  and  numerically,  and  in 
the  case  of  vortex  sheets  with  no  surface  tension  certain  solutions  are  given  in  closed  form. 


'This  research  was  supported  by  the  National  Aeronautics  and  Space  Administration  under  NASA  Con¬ 
tract  Nos.  NAS1-I9480  and  NASl-1860.5  while  the  author  was  in  residence  at  the  Institute  for  Computer 
Applications  in  Science  and  Engineering  (IC-'ASE),  NASA  Langley  Research  Center,  Hampton,  VA  2.3681. 

^This  research  was  supported  in  part  by  Fondo  Nacional  de  Desarrollo  Caentifico  y  Tecnologico  (FONDE- 
CYT)  under  grant  No.  246/91  and  Universidad  Tecnica  Federico  Santa  Maria,  Valparaiso,  t^hile. 


Ill 


1  Introduction 

Rayleigh  [1]  first  studied  the  break-up  of  liquid  jets  in  air.  The  underlying  instability  mech¬ 
anism  is  a  consequence  of  surface  tension,  and  Rayleigh  predicted  that  all  waves  longer  than 
the  jet  diameter  are  linearly  unstable  and  used  the  result  to  conjecture  the  break-down  of 
the  jet  into  drops.  The  physics  of  the  instability  can  be  understood  in  purely  energetic 
terms  since  for  a  given  stationary  volume  of  fluid  the  surface  energy  is  a  minimum  when 
the  shape  is  spherical.  A  liquid  thread,  therefore,  will  tend  to  break  up  into  drops.  An 
equivalent  and  more  quantitative  feature  of  the  instability  is  found  also  by  consideration  of 
the  Young-Laplace  equation  which  relates  the  pressure  jump  across  the  jet  interface  to  its 
curvature.  In  cylindrical  geometries  we  have  two  piiucipal  radii  of  curvature,  one  in  a  plane 
containing  the  axis  of  the  jet  and  one  in  a  plane  perpendicular  to  the  axis.  If  we  suppose 
that  a  perturbation  causes  the  jet  radius  to  decrease  slightly,  then  the  radius  of  curvature 
in  the  cross-sectional  plane  decreases  and  the  pressure  in  the  vicinity  of  the  perturbation 
increases.  This  means  that  fluid  will  tend  to  leave  this  high  pressure  region  to  flow  to  lower 
pressure  positions,  which  by  mass  conservation  decreases  the  jet  radius  further  thus  enhanc¬ 
ing  the  interfacial  instability.  This  is  a  physical  description  of  the  mechanism  of  pinching 
or  break-up.  Clearly  the  phenomenon  is  nonlinear  and  even  though  linear  theory  does  well 
in  predicting  dominant  wavelengths  in  natural  transition,  nonlinear  analyses  are  required  to 
follow  the  jet  to  break-up  for  a  range  of  dominant  length  scales  whether  self-sustained  or 
externally  imposed  through  initial  conditions  for  example. 

The  physical  picture  sketched  above  is  fundamental  and  appears  in  both  inviscid  as  well 
as  viscous  flows.  In  this  work  we  concentrate  on  inviscid  irrotational  flows  of  liquid  jets.  A 
review  article  has  been  written  by  Bogy  [2].  Of  relevance  is  also  the  work  of  Chandrasekhar 

[3]  and  some  experimental  pictures  of  jet  break-up  can  be  found  in  Donnelly  and  Glaberson 

[4]  (see  also  Drazin  and  Reid  [5]).  A  theoretical  study  of  the  break-up  of  an  inviscid  incom¬ 
pressible  liquid  jet  in  air  with  surface  tension  has  been  carried  out  by  Ting  and  Keller  [6] 
(referred  to  as  [TK]).  The  authors  make  a  slender  jet  approximation  where  the  radial  length 
scale  is  taken  to  be  much  smaller  than  the  axial  length  scale  so  that  their  ratio  forms  a  small 
parameter  which  is  used  in  an  asymptotic  theory  to  develop  a  system  of  one-dimensional 
coupled  PDE’s.  Similar  ideas  have  also  been  used  by  Papageorgiou  and  Smith  [7]  in  the 
study  of  the  nonlinear  evolution  of  plane  wakes.  The  evolution  equations  of  [TK]  are  simpler 
than  the  original  system  making  them  amenable  to  an  analysis  that  produces  families  of 
pinching  solutions  which  are  physically  relevant.  Self-similar  solutions  in  turn  transform  the 
PDE’s  into  ode’s  and  numerical  or  exact  solutions  are  given. 

The  problem  studied  here  is  an  extension  of  [TK].  In  [TK]  the  cylindrical  jet  flows  in  air. 
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Within  the  framework  of  the  Euler  equations  the  undisturbed  velocity  field  is  a  constant  axial 
velocity  which  can  be  removed  from  the  problem  by  a  suitable  Galilean  transformation.  The 
mechanism  of  three-dimensional  Kelvin-Helmoltz  instability  is  absent  and  surface  tension 
dominates.  Our  interest  lies  in  the  competing  effects  of  capillary  instability  and  Kelvin- 
Helmholtz  instability  which  is  present  if  the  surrounding  (without  loss  of  generality  ambient) 
fluid  is  of  non-zero  density.  Recently  there  has  been  a  lot  of  interest  in  the  formation  of 
singularities  in  cylindrical  vortex  sheets.  Linear  stability  studies  include  the  work  of  Rotunno 
[8]  and  more  recently  singularity  formation  studies  have  been  undertaken  by  Caflisch  [9], 
Caflisch  and  Li  [10],  Caflisch,  Li  and  Shelley  [11],  Grauer  and  Sideris  [12],  who  do  not 
include  surface  tension  but  include  a  background  swirl  in  the  azimuthal  direction.  The  flows 
in  [9]-[ll]  do  not  possess  a  jump  in  axial  velocity  at  the  sheet  surface  and  are  shown  in 
[8]  and  [11]  to  be  linearly  stable  to  short  waves  in  the  case  of  axisymmetric  disturbances. 
The  problem  studied  here  omits  swirl  and  includes  axial  velocity  jumps  which  contribute  a 
Kelvin-Helnilioltz  instability  for  axisymmetric  disturbances.  Axisymmetric  flows  are  more 
likely  to  form  pinching  singularities  since  the  axisymmetry  prevents  the  core  from  flattening 
out.  Even  though  we  do  not  include  swirl  in  our  study  we  do  enforce  axisymmetry,  and  our 
slender  jet  equations  predict  self-similar  pinching  of  the  jet  even  in  the  absence  of  surface 
tension.  Such  solutions  may  be  suitable  candidates  for  initializing  full-scale  computations. 

The  rest  of  the  paper  is  organized  as  follows.  Section  2  derives  the  slender  jet  equations. 
For  completeness  we  derive  a  system  for  fluids  of  different  densities  with  the  axis  of  the  jet 
being  vertical.  There,  we  also  describe  some  useful  transformations  that  reduce  the  num¬ 
ber  of  unknowns  in  the  evolution  equations  from  three  to  two  which  are  also  more  relevant 
physically.  Section  3  discusses  properties  of  the  equations  including  conserved  quantities, 
equation  type  and  linear  stability  characteristics.  According  to  linear  theory,  stability  win¬ 
dows  exist  for  certain  values  of  density  ratios  and  certain  undisturbed  jet  to  pipe  radius 
ratios.  Equal  density  flows  are  always  linearly  unstable.  In  Section  4  we  consider  the  possi¬ 
bility  of  finite-time  singularities  and  derive  the  self-similar  governing  equations.  In  Section  5 
we  solve  these  equations  for  vortex  sheets  in  the  absence  of  surface  tension  and  give  families 
of  exact  solutions,  and  in  Section  6  we  address  the  problem  including  surface  tension.  Section 
7  is  devoted  to  discussion. 

2  Derivation  of  the  evolution  equations 

Consider  the  irrotational  flow  of  an  axisymmetric  jet  of  density  p\  and  potential  4>]{t,r,z) 
whose  bounding  surface  is  S{t,z).  The  flow  outside  is  also  irrotational  and  is  taken  to 
have  potential  (l>2{t,r,z)  and  density  p^-  The  flow  takes  place  in  a  vertical  arrangement 
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with  gravity  acting  and  there  is  an  outer  bounding  solid  surface  at  r  =  R-i  (a  cylindrical 
polar  coordinate  system  (r,  6,  z)  is  used  throughout  with  the  corresponding  velocity  field 
u  =  (u,  u,io)).  If  in  addition  surface  tension  acts  and  the  tension  coefficient  is  constant,  cr 
say,  the  following  piecewise  constant  solution  of  the  Euler  equations  is  available: 


Ui  =  (0,0,  lio)  0  <r  <  Ri, 


u-i  =  (0,0,0)  R%  <r  <  Ri, 

where  Wo,  R\  are  constants.  This  solution  represents  a  perfectly  cylindrical  jet  of  constant 
axial  velocity  flowing  into  a  quiescent  medium.  When  the  jet  interface  deforms  we  need  to 
solve  the  problem  in  the  regions  Vi  :  0  <  r  <  S{t,z)  and  T>2  :  S{t,  z)  <  r  <  Ri.  There  are 
two  separate  mechanisms  of  instability  present:  (i)  Kelvin-Helmholtz  instability  of  a  cylin¬ 
drical  vortex  sheet,  (ii)  the  destabilization  of  long  waves  by  the  circumferential  component 
of  the  radius  of  curvature  if  surface  tension  acts.  We  describe  the  nonlinear  evolution  by 
assuming  that  the  axial  length  scale  of  interfacial  perturbations  is  large  compared  to  the 
radial  length  scale.  We  note  that  such  long  length  scales  do  not  arise  directly  from  linear 
stability  analyses  of  jets  (the  characteristic  length  is  of  the  order  of  the  jet  radius  for  capillary 
dominated  flows),  but  enter  the  nonlinear  evolution  through  initial  interfacial  conditions  of 
the  form  /{(z)  with  e  small.  The  following  non-dimensionalizations  are  appropriate: 


--_L  2  -_L 

'  Ri  ^  d  ‘  R^  ^  dWo 


<^1,2 


(2.1) 


Utilization  of  the  scalings  (2.1)  into  the  governing  equations  (omitting  the  bars)  gives: 
Laplace  equations  in  V^^i 

<f>rr  H - <t>r  +  =  0  (2-2) 

V 

Kinematic  conditions  on  r  =  S 


<f>\r  —  (^{St  -f  <j>\z^z) 
hr  =  e'(S',  -b  h.S,) 


Bernoulli  equation  on  r  =  S 


hi  +  +  4>\z)  -  Pihi  + 


=  -w. 


(1  -b 


s 


(1  -b  f^5?)V2  (1  +  fKSi 


(2.3a) 

(2.36) 


(2.4) 
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Boundary  condition  r  =  a 


(t>2T{r  =  a)  =  0. 


(2.5) 


The  constant  terms  in  the  Bernoulli  equation  (2.4)  are  chosen  to  satisfy  the  unperturbed  flow. 
The  term  on  the  right  hand  side  is  the  increment  induced  by  surface  tension  on  the  pressure 
acting  across  the  interface.  In  (2.2)-(2.5)  above,  the  following  nondimensional  groups  appear 


P2 


1 


--p  ^-^(“  =  TT7' 


w  _  ^ 

^  R,p,Wi  “  R, 


These  parameters  are  the  ratio  of  radial  to  axial  length  scale,  the  density  ratio  between  outer 
and  inner  fluid,  the  Weber  number  (ratio  between  capillary  pressure  and  dynamic  pressure), 
and  the  ratio  between  the  cylinder  radius  and  the  unperturbed  jet  radius.  Equations  (2.2)- 
(2.5)  are  exact  and  aie  a  result  of  a  change  of  variables.  For  slender  jets,  however,  e  is  a 
small  parameter  and  the  following  asymptotic  expansions  are  appropriate 


4>i;2  —  +  .  • . ,  (2.6a) 

5  =  So  +  +  . . .  .  (2.66) 

Substitution  of  (2.6a,b)  into  (2.2)  gives  the  following  solutions  that  satisfy  boundedness  at 
the  axis  (r  =  0)  and  the  boundary  condition  at  r  =  o 


,  <l>P  =  B2{t,z),  (2.7) 

4>\  ^  =  —^B^zz  +  C]{i,z)  , 

^1)  =  -jB^z.  +  jB-izzlnir)  +  C2{t,z).  (2.8) 

The  functions  5i,2  and  (7j,2  are  to  be  determined  from  the  kinematic  and  Bernoulli  equations 
(2.3)  and  (2.4).  The  solutions  (2.7)  for  the  leading  order  potentials  in  Pi, 2  automatically 
satisfy  (2.3a,b)  and  (2.4)  to  leading  order.  Next,  at  O(f^)  and  0(1)  respectively,  equations 
(2.3a,b)  and  (2.4)  along  with  the  solutions  (2.7),  (2.8)  give 


~'^SqBizz  —  Sot  +  BtzSozi 
-UoB-izz  +  ^B2zz  =  Sot  +  B^zSoz, 

Z  Z*70 

1  1  1 

ei.  +  t^b],  -  p{B„  +  -Bl)  =  —^  +  W,- -Wl 


(2.9) 

(2.10) 

(2.11) 


Equations  (2.9)-(2.11)  as  given  above  constitute  a  nonlinear  system  of  three  evolutions  equa¬ 
tions  for  the  three  unknowns  .Sq,  B\,  Bi.  Next  we  show  how  they  can  be  reduced  to  a  .set  of 
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two  coupled  equations  for  the  shape  of  the  jet  and  the  axial  velocity  jump  across  it.  This 
reduction  is  useful  both  mathematically  and  physically  in  discussing  the  pinching  process 
after  a  finite  time  (see  later). 

If  we  subtract  (2.10)  from  (2.9)  to  eliminate  .S'ot,  the  result  can  be  readily  integrated  once 
with  respect  to  z  and  enables  expression  of  in  terms  of  by  the  formula 

=  (2.12) 


The  unknown  function  of  t  that  arises  from  the  integration  can  be  taken  to  be  identically 
equal  to  zero  by  choice  of  appropriate  initial  conditions  (see  Papageorgiou  and  Smith  [7]  also 
for  a  similar  theory  for  two-dimensional  wakes).  Next  we  define  new  variables  W  and  V 
given  by 

W  =  B,-B2  ,  V=B,  +  B2  ,  (2.13) 

which  transform  (2.12)  into 

.  (2.14) 

The  first  evolution  equation  is  found  by  addition  of  the  kinematic  conditions  (2.9)  and  (2.10), 
and  elimination  of  V  while  the  second  equation  follows  from  (2.11)  again  after  elimination 
of  V  using  (2.14).  It  is  convenient  to  differentiate  the  resulting  Bernoulli  equation  with 

respect  to  z  for  two  reasons;  (i)  The  constant  (which  can  also  be  a  function  of  time  in  other 

configurations)  associated  with  the  Bernoulli  integral  disappears,  (ii)  the  resulting  equations 
involve  the  functions  5o  and  Wz.  The  latter  is  seen  to  be  the  jump  in  axial  velocity  across 
the  deforming  vortex  sheet  and  is  therefore  proportional  to  the  sheet  strength.  Introducing 
the  variables 

r/  =  %  ,  X  =  Wz  , 

gives  the  following  coupled  system  of  evolution  equations; 

m  +  ivx  -  v^x)z  =  0’ 

+  «(^  -  »?)]xt  +  otx{vx  -  v^x)^ 

+(l|  +  f  (i  i(i  -  ,)1X=)^  =  ;  P  ”) 

Equation  (2.17)  above,  simplifies  considerably  for  the  case  of  equal  densities  a  =  0;  equation 
(2.16)  is  unaffected  by  density  differences  and  (2.17)  for  ft  =  0  becomes 


(2.15) 


(2.16) 


(2.18) 


Equations  (2.16)  and  (2.17)  (or  (2.18))  are  to  be  solved  subject  to  prescribed  conditions  at 
<  =  0.  To  fix  matters,  we  will  assume  that  the  flow  is  periodic  in  of  period  27r  (a  flow  of 
any  other  period  can  be  scaled  to  this  canqnical  case  by  rescaling  time). 


3  Properties  of  the  evolution  equations 

3.1  Conserved  quantities 

The  following  are  conserved  quantities  of  the  evolution 

fiir  J  1 

^  x{-2  +  2^  -  ,  ^2  =  y^  r]dz.  (3.1a, 6) 

The  first  conserved  quantity  is  a  direct  consequence  of  the  Bernoulli  equation  (2.11).  The 
total  energy  of  the  system  is  also  a  conserved  quantity.  The  calculations  that  lead  to  (3.1c) 
below  are  included  in  Appendix  A. 

f  /\  1  W  .  \ 

^3  =  j^  [jp{i J  (3-lc) 


3.2  Equation  type 

Equations  (2.16)  and  (2.17)  can  be  written  in  the  form 

where  A  is  a  2  x  2  matrix  with  entries  which  are  nonlinear  functions  of  tj  and  x-  It  Is  easy 
to  show  that  the  eigenvalues  of  A,  A  say,  satisfy 

(x-2,x-a)((1±^i:M_a)  = 


^0 _ ^ 

A 


( 


l)x^ 


2a(l  +  p) 


where 

*  I  A 
A=  2  +  0(2--)). 

The  type  of  the  evolution  equations  depends  on  the  roots  of  the  polynomial  above  !n  what 
follows  we  show  that  the  eigenvalues  are  complex,  in  almost  all  cases,  and  so  the  system  is 
elliptic  in  general.  To  do  this  it  is  enough  to  show  that  the  discriminant  of  the  quadratic 
equation  is  negative.  After  a  series  of  algebraic  manipulations  it  is  shown  that  the  eigenvalues 
are 

where 

jj2  ^  4x^?/(l  -  ?7)(a^  -  1)  _  4M47?(1  - 

(1  d-  a  —  2aTjy  a(l  +  p)(l  +  «  —  2arf) 
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Since  0  <  t;  <  1  and  <  1  we  can  conclude  that  the  eigenvalues  are  complex  since  <  0. 
The  reason  for  this  is  that  both  terms  in  the  expression  for  above  are  negative;  this 
conclusion  uses  the  fact  that  2A  =  (1  +  a  —  2aj;)  >  0  for  all  relevant  values  of  a  and  r}. 

It  is  concluded,  therefore,  that  the  system  is  elliptic  for  almost  all  values  of  o  and  W^. 
The  exceptional  case  has  a  =  1,  =  0  which  represents  a  jet  in  air  with  zero  surface  tension 

yielding  real  but  repeated  eigenvalues  (further  discussion  of  this  case  is  given  in  Section  7). 

The  linear  stability  spectrum  follows  from  a  linearization  of  the  system  (3.2)  about  the 
unperturbed  states.  If  we  write 

it  follows  that  c  is  given  by  ±{D^y^^  where  Z)  is  as  given  by  (3.3)  and  evaluated  at  {x^v)  = 
(1,^).  As  shown  above,  <  0  and  so  instability  is  obtained  with  growth  rates  of  0{k), 
as  in  Kelvin-Helmholtz  instability  of  a  flat  ''ortex  sheet  for  instance.  Our  aim  here  is  to 
follow  this  instability  into  a  fully  nonlinear  regime  with  the  cylindrical  jet  interface  growing 
to  amplitudes  comparable  to  the  unperturbed  radius. 


4  Finite-time  singularities:  Jet  pinching 


In  this  section  we  consider  analytically  possible  nonlinear  terminal  states  of  the  evolution 
equations.  It  is  found  that  density  differences  are  unimportant  in  the  leading  order  balances 
between  unsteady  and  nonlinear  terms  as  a  singularity  is  encountered.  It  is  sufficient  to 
consider  equal  densities  p  =  1  (a  =  0).  The  case  p  ^  I  also  leads  to  equations  (4.4a,b) 
below,  with  Wg  replaced  by  'IWg.  For  equations  (2.16)  and  (2.18),  therefore,  we  postulate 
that  the  solutions  encounter  singularities  after  a  finite  time,  i,  say,  at  axial  positions  Zg  and 
seek  to  describe  the  solutions  near  such  points.  We  make  the  following  ansatz 


77(2,  t)  =  (t,  - 1)“/(0  x(^,0  =  {t.  -  tVg{i) 

f  =  tr  :■) .  (4,1) 

Physically  relevant  solutions  have  a  and  (i  positive  in  order  to  ensure  that  the  inLfifai.e 
position  does  not  become  unbounded  at  the  critical  time,  and  that  the  dominant  features 
are  due  to  nonlinear  focussing.  Under  the  transformations  (4.1)  time  and  space  derivatives 
become 


±  9  (H  ^  1  9 

dt^  {tg-t)d^  dz~"  ^  ’ 

Substitution  of  (4.1),  (4.2)  in  (2.16)  and  (2.18)  along  with  the  usual  order-of-magnitude 
arguments  as  t  from  below  gives  the  following  expressions  relating  the  constants; 


a  =  4  —  40  7  =  /Z  —  1 , 


(4.3) 
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and  for  a  to  remain  positive  we  require 

0<.j3<L 

Note  that  this  choict.  ol  implies  a  negative  7  and  so  the  axial  velocity  jump  blows  up  as 
the  jet  radius  shrinks  to  zero. 

The  conserved  integrals  (3.1a-c)  become  to  leading  order  in  (t,  —  i): 

/OO  /•OO 

gd(  1-2  =  iu  - /  fd^ 
h  =  («.  -  ()"■"  /“(t/Z  +  —/’'’)  rff-  ('•■'I) 

J-00  4  a 

Equations  (4.4)  above  represent  the  local  contributions  of  the  focussing  region  to  the  three 
conserved  quantities.  If  the  exponents  of  {t,  ~  t)  are  positive  the  contributions  are  zero 
in  the  limit  and  no  integral  constraints  need  to  be  placed  on  possible  solutions  /(^)  and 
5r(^),  whereeis  if  any  exponent  is  negative  or  zero  the  corresponding  integral  it  multiplies 
must  vanish.  Using  (4.3a)  we  see  that  the  I2  and  I3  contributions  are  zero  in  the  limit. 
Consideration  of  the  /j  contribution,  however,  shows  that  the  following  constraint  is  required: 

0<\/2  rgd^  =  0.  (4.5) 

J  — OO 

The  governing  equations  transform  to  the  following  one-parameter  system  of  nonlinear 

ode’s: 

-4{l-f3)f  +  /)a'  +  {fgy  =  0,  (4.6n) 

(1  -  Hh  +  /%'  +  (4.61.) 

Z  a 

In  general  (4.6a,b)  must  be  solved  numerically  for  different  values  of  fi  and  with  conditions 
at  ^  =  0  (i.e.  z  =  z,)  specified  for  example. 

Before  the  presentation  of  solutions  to  (4.6a,b),  we  note  that  the  same  self-similar  equa¬ 
tions  are  obtained  for  the  case  of  a  jet  in  air.  Further  more,  the  constraints  on  imposed 
by  the  conserved  integrals  of  the  system  are  preserved  and  so  the  sohitions  that  follow  ap¬ 
ply  equally  well  to  jets  in  air.  For  completeness  the  analogous  derivation  of  a  jet  in  air  is 
contained  in  Appendix  B. 

5  Collapse  of  cylindrical  vortex  sheets  with  zero  sur¬ 
face  tension. 

The  equations  governing  the  self-similar  break-up  structure  become: 

-4( l-/^)/  +  /y^/'-h  (/(/)' =  0,  (5.1«) 
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(l-/3)5  +  /i^(7'  +  -(/)'  =  0.  (5.16) 

The  equation  for  g  is  nonlinear  but  is  decoupled  from  that  for  /  which  is  linear.  Once  g  is 
determined,  the  solution  for  /  readily  follows; 


UH  +  gy 


(5.2) 


The  main  task,  therefore,  is  in  the  determination  of  g  and  classes  of  exact  solutions  are 
constructed  next. 

An  exact  solution  is  readily  available  for  the  critical  case  (5.1b)  can  be  integrated 

once  to  yield  a  quadratic  algebraic  equation  for  <7(^)  which  yields  the  solutions 


2<7±  =  -^  ±  +  4^^, 


(5.3) 


where  =  </(0)-  The  solution  g+  pertains  to  positive  velocity  jumps  while  </_  describes 
flows  with  negative  jumps.  The  corresponding  solutions  for  /  are 


A  = 


fo  {i+\/e+‘id)‘ 
\6gi  (e+4s2)’'“ 


/-  =  M/osiu  +  ]/e+‘igSrH(‘ + 

with  /o  =  /(O).  The  behavior  for  large  ^  is 

fo 


(5.4) 


U 


*  ^ 

4  s  i7+ 


16</o 

U'^2gl\i\-^  g^^ 

/_~64^oVor"  g- 
f- 


+00, 


—1^1“* 

9o 


fir-  ~  - 


kl  ^  -oo- 

^  +0O, 

^  —*  — CX). 


gl 

1^1 


Even  though  the  solutions  f±{Q  are  admissible  the  corresponding  ^±(0  i^ot,  due  to 
the  conservation  constraint  (4.5)  requiring  that  the  integral  vanish.  This  is  satisfied  by 
<7(0  =  fir- (0  for  4  <0,  andi?(0  =  <7+(0for^  >  0,  but  such  solutions  are  not  acceptable  since 
g{^)  and  /(O  are  now  discontinuous  at  ^  =  0,  in  violation  of  our  slender  jet  approximation. 
We  reiterate  here,  that  any  solutions  /(O?  g{0  bounded  derivatives  are  acceptable 

singular  states  of  the  slender  jet  equations;  unbounded  derivatives  herald  the  introduction 
of  shorter  scales  over  which  a  new  acceptable  solution  may  emerge.  Such  analyses  are  not 
pursued  further  here  since  we  are  interested  in  the  construction  of  pinching  solutions. 

For  any  value  of 

g(i)  =  -(  HO  = 
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are  exact  solutions  with  A  some  constant.  The  solution  for  g  is  an  odd  function  of  ^  ami  is 
consistent  with  the  integral  constraint  (4.5).  For  /(^)  to  he  acceptable  we  require  bounded¬ 
ness  of  the  first  derivative  at  =  0;  this  is  not  possible  for  ft  <  \,  and  this  family  of  solutions 
is  therefore  dropped. 

For  ft  =  2/3  another  exact  solution  is  found:  this  is  most  easily  seen  by  multiplication  of 
(5.1b)  by  g  and  integration  with  respect  to  ^  which  gives 

+  '-Kg‘  + 1  =  ao  =  5(0), 

from  which  it  follows  that  for  ft  =  2/3,  g  satisfies  the  cubic  polynomial 

+  =  (5.5) 

If  <7o  =  0  the  nontrivial  solution  is  </  =  — ^  with  corresponding  /  =  -4^“^  which  is  singular 
at  the  origin  and  is  discarded  (see  above).  Equation  (5.5)  can  have  (i)  one  real  root  and  a 
pair  of  complex  conjugate  roots,  (ii)  three  real  roots  of  which  at  least  two  are  equal,  and  (iii) 
three  distinct  real  roots,  depending  on  the  conditions 


(i)  Q  >  0  (ii)  Q  =  0  (Hi)  Q  <0, 


Q  =  J<7o(<7o  - 

respectively.  The  position  ^  ~  =  '^9o  yields  (5  =  0  and  is  the  point  above  which  three 
distinct  solutions  exist.  At  ^  =  there  is  a  root  g  =  (7o/4’'^^  and  a  pair  of  equal  roots 
g  =  —2go(4^^^.  A  local  analysis  shows  that  <7  ~  (^  —  1)’''^^  in  the  vicinity  of  the  double  root 
indicating  that  the  derivative  is  infinite  at  The  solutions  of  (5.5)  for  gQ  =  4/27  are  shown 
in  Figure  1;  the  dashed  line  represents  g  =  —^  which  is  the  asymptotic  state  of  one  of  the 
solutions  for  large  Within  the  context  of  the  long  wave  theory,  the  roots  with  unbounded 
derivatives  are  discarded  and  so  a  unique  acceptable  solution  exists  for  ft  =  2/3. 

The  type  of  solution  given  above  for  ft  =  2/3  can  be  generalized  to  arbitrary  values  of  ft. 
Under  the  transformation 

9iO  =  (5.6) 


equation  (5.1b)  becomes 

G’^(l+G)'-^  =  y,  (5.7) 

where  is  a  constant  of  integration.  Equation  (5.7)  yields  G  implicitly  in  terms  of  We 
can  rewrite  (5.7)  in  terms  of  g  using  (5.6)  to  give 


+  5')’  ^  —  9oy 


(5.8) 
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where  go  =  <7(0).  Equation  (5.8)  is  a  generalization  of  (5.5)  to  arbitrary  values  of  (S.  In 
particular  if  /?  is  a  rational  number  with  ft  —  pjq  with  p,q  integers  such  that  Ij'l  <  ft  <  1, 
(5.8)  becomes 

9’’«  +  jr"  =  sS-  (5!>) 

The  values  p  =  2,  q  =  'i  yield  equation  (5.5).  We  are  interested  in  generating  real  solutions 
of  (5.9)  which  are  acceptable  similarity  forms  of  the  breakdown.  To  do  this  we  begin  at 
=  0  where  there  is  one  real  root  and  <7  —  1  complex  roots  if  q  is  odd,  or,  two  equal  and 
opposite  real  roots  and  q  —  2  complex  roots  if  q  is  even,  and  follow  the  appropriate  real  root 
as  a  function  of  (  to  obtain  g{^)-  This  can  be  done  for  negative  ^  also;  alternatively  the 
differential  equations  can  be  solved  directly.  This  was  also  done  as  a  check. 

The  behavior  of  solutions  to  (5.1a,b)  for  large  ^  is  found  to  be 

(^00.  (5.10a) 

The  constants  of  proportionality  can  be  determined  numerically.  For  ^  — >  —00  the  equations 
are  analyzed  for  large  positive  ^  to  give  the  asymptotic  forms 

/~(-0"’  5  =  (-0  +  ---  ^ (O-106) 

We  note  that  the  constant  of  proportionality  for  the  decay  rate  of  /  is  determinable  numeri¬ 
cally,  whereas  the  constant  for  the  leading  order  asymptotic  behavior  for  g  is  equal  to  unity. 
These  results  are  confirmed  by  our  numerical  solutions  described  next. 

Numerical  solutions  can  be  generated  for  any  ft  in  the  designated  range.  Here  we  present 

a  set  of  solutions  for  ft  =  2/3,4/5,9/10.  The  results  are  presented  in  Figures  2,  3  for  /(O 

and  g{(^)  respectively.  The  initial  conditions  for  all  solutions  are  /(O)  =  0.5,  g{0)  =  1.0. 
The  asymptotic  results  (5.10a,b)  are  confirmed  by  these  solutions.  The  behavior  of  g{()  as 
^  — >  —00  predicted  above  is  clearly  seen  in  Figure  8  where  the  behavior  is  independent  of 
ft  and  the  curves  collapse  to  a  line  of  unit  slope.  Of  particular  interest  are  the  solutions  for 
ft  =  4/5.  According  to  (5.10a),  /(^)  grows  linearly  for  large  positive  ^  when  ft  =  4/5.  For 
ft  <  4/5  the  f{(ft  is  concave  for  large  ^  and  for  ft  >  4/5  it  is  convex.  As  explained  below, 
viewed  from  outside  the  similarity  region,  the  solution  for  ft  =  4/5  represents  a  conical  jet 
(the  radius  is  a  linear  function  of  axial  position)  while  the  solutions  for  ft  <  4/5  and  ft  >  4/5 
repre.sent  concave  and  convex  jets  respectively  at  the  singular  time.  The  convex  jet  solution 
which  has  ft  =  7/8  is  shown  in  Figure  4  for  ^  >  0.  The  initial  conditions  are  /(O)  =  0.3, 
^(0)  =  0.6.  The  a.symptotic  forms  (5.10a)  are  confirmed  numerically  in  Figure  5  by  use  of 
a  logarithmic  plot  and  a  least  squares  fit  of  the  resulting  data  to  obtain  an  estimate  for 
the  slope.  As  expected  the  errors  between  the  exact  and  the  computed  value  of  the  slope 
decrease  slowly  as  the  upper  limit  of  integration  is  increased. 
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Since  the  solution  for  /  becomes  unbounded  for  large  positive  ^  (as  does  that  for  g  at 
large  negative  ^),  it  is  essential  to  check  that  these  asymptotic  values  are  consistent  with 
the  ansatz  (4.1);  such  analysis  also  provides  a  qualitative  description  of  the  break-up  from 
outside  the  similarity  region.  If  we  substitute  the  asymptotic  forms  (5.10a)  into  (4.1)  and 
express  ^  in  terms  of  the  physical  spatial  variable  Z-Zs{>  0)  we  find  that  when  viewed  from 
outside  the  self-similar  region,  the  break-up  takes  the  spatial  characteristics 

Tj  ^{Z-  Zs)  X  ~  (-  -  -,)  . 


For  z  <  Zg  the  asymptotic  forms  (5.10b)  are  used  to  give 

{tg  -  (.=:,  -  2 

[Zg  -  2)  (tg  -  f. 

These  expressions  indicate  the  geometry  of  the  jet  at  break-up;  the  shape  is  squeezed  to  a 
thread  as  ^  >  —00  with  the  radius  tending  to  zero  as  the  singular  time  is  approached  (this 

cai]  be  seen  from  the  fact  that  1  <  4  —  3/i  <  5/2)  on  one  side,  while  it  has  a  concave  or 
convex  growth  on  the  right  depending  on  the  value  of  /i.  The  corresponding  velocity  jumps 
tend  to  infinity  and  zero  respectively  at  the  algebraic  rates  given  above.  For  example,  the 
exact  solution  given  above  for  /i  =  2/3  gives  a  parabolic  shape  at  break-up. 


6  Collapse  of  cylindrical  vortex  sheets  including  sur¬ 
face  tension 


The  governing  equations  are  (4.6a,b).  The  parameter  W^/a  in  (4.6b)  ran  be  re-scaled  to 
unity  by  the  transformation  /  iW^fa^)f,  yielding  the  normalized  form 

-4(  I -/^)/-b/^a' +  {/</)' =  0,  (6.1a) 

(1  -  /%  +  «/+)(/)'  =  -(/-■'")'.  ilM) 

Before  presentation  of  numerical  solutions  it  is  useful  to  examine  the  large  |^|  behavior  of 
the  solutions.  For  large  positive  ^  we  look  for  solutions  of  the  form 

/  ~  AC  9  ~  Hr"  rn,  n  >  0.  (6.2) 


Using  these  expressions  in  (6.1a)  implies 


m  = 


4(1  -/j) 
/i 


12 


The  correct  value  of  n  which  matches  the  numerical  solutions  is  found  by  a  balance  of  the 
linear  terms  in  (6.1b).  To  illustrate  this  we  write  the  asymptotic  form  of  (6.1b)  according  to 
the  scalings  (6.2), 

(1  -  ll)Br‘  -  linBC"  -  e'nr'”-'  + 

with  part  of  the  higher  order  correction  terms  shown  explicitly.  If  we  choose  to  match  the 
linear  terms  in  B  on  the  left  with  the  surface  tension  nonlinearity  on  the  right,  a  value  for 
71  is  found  which  makes  the  constant  B  negative  in  contrast  to  numerical  solutions.  The 
correct  behavior  is  determined  by  matching  the  linear  terms  which  gives 

_  1  — 

”  ^  T  " 

With  this  value  of  n  the  linear  terms  (which  are  the  largest  terms  in  the  equation)  match 
automatically  with  B  arbitrary.  Similar  arguments  apply  for  ^  — +  — oo  (the  analysis  is 
almost  identical  if  we  write  ^  and  analyze  the  equations  as  +  oo).  To  summarize, 

therefore,  we  find 

^-^oo  (6.:3a) 

(6.36) 

Note  that  the  solutions  for  large  positive  ^  have  the  same  asymptotic  behavior  as  those  corre¬ 
sponding  to  zero  surface  tension.  As  ^  >  — oo,  however,  the  behavior  is  altered  significantly. 

The  consequences  of  this  are  discussed  after  the  presentation  of  numerical  solutions. 

Equations  (6.1a,b)  were  solved  numerically  with  particular  interest  in  the  verification  of 
the  asymptotic  structures  (6.3a,b).  The  initial  conditions  used  were  /(O)  =  1/2,  g{0)  = 
1  with  an  integration  range  as  large  as  (—100,100).  Results  are  presented  for  the  three 
representative  values  fi  =  2/3,4/5,9/10.  Figure  6  shows  collectively  the  solutions  for  g{^). 
The  maximum/minimum  in  g  decreases/increases  as  0  is  increased.  A  notable  feature  is 
that  there  are  two  distinct  regions,  ^  <  >  4o  for  some  value  where  the  axial  velocity 

jump  is  negative  and  positive  respectively.  Changing  the  initial  condition  so  that  (/(O)  <  0 
for  example,  indicates  that  this  separation  is  preserved.  Results  for  the  corresponding  / 
solutions  are  given  in  Figures  7a,b,c  separately  due  to  the  large  scale  differences  involved. 
For  /3  <  4/5  the  computed  shape  is  concave,  and  the  solution  for  /i  =  2/3  has  a  minimum  at 
^  =  —0.3  approximately;  the  solution  corresponding  to  /i  =  4/5  is  a  critical  one  with  a  linear 
growth  at  ±oo  (see  Section  5  also),  with  a  minimum  radius  at  ^  =  —0.38  approximately; 
solutions  having  0  >  4/5  are  convex  in  regions  sufficiently  far  away  from  the  origin  and 
the  computed  solution  for  0  =  9/10  has  a  minimum  at  ^  =  —0.48  approximately.  For  the 
solutions  presented  here,  the  position  of  minimum  radius  (i.e.  the  value  of  ^  which  makes 
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/(^)  a  global  minimum)  is  larger  than  the  positions  where  the  jump  in  axial  velocity  (which 
proportional  to  g{^))  is  zero. 

Next  we  use  the  numerical  solutions  just  presented  to  check  the  2isymptotic  forms  of  the 
solutions.  This  is  done  using  logarithmic  plots  and  a  least  squares  fitting  of  the  resulting 
data.  For  large  positive  (  results  are  shown  in  Figures  8a-c  for  /  and  g  for  the  three 
representative  values  of  (j.  The  agreement  is  very  good  and  is  found  to  improve  as  the  value 
of  computational  “infinity”  is  increased,  as  expected.  For  large  negative  ^  analogous  results 
with  similar  accuracy  are  found  and  need  not  be  presented  here. 

As  in  Section  5,  the  asymptotic  form  of  the  solutions  far  from  the  origin  enables  the 
qualitative  description  of  the  shape  at  break-down.  Considering  the  shape  first,  there  appear 
to  be  three  canonical  forms  at  break-down  corresponding  to  the  cases  (i)  /5  <  4/5,  (ii) 
0  =  4/5,  and  (iii)  >  4/5.  For  case  (i)  a  concave  shape  is  found  at  infinity  and  the  jet 
would  tend  to  pinch  forming  drops;  case  (ii)  gives  linear  shapes  at  infinity  and  at  break¬ 
down  the  jet  surface  would  look  like  two  right  circular  cones  with  a  common  vertex;  case 
(iii)  yields  jets  with  shapes  which  are  locally  convex  in  the  axial  direction.  Photographs  of 
experiments  in  air  (see  [4],  [5])  indicate  behavior  which  fits  in  with  cases  (i)  and  (iii).  An 
additional  fact  which  is  evident  from  the  experiments  and  makes  long  wave  theories  like  the 
present  one  (see  also  [TK])  useful  candidates  for  break-down  prediction,  is  the  validity  of  the 
slender  jet  assumption  and  in  particular  the  absence  of  gradient  steepening  in  the  shapes  at 
break-down. 

7  Conclusions 

We  have  developed  a  nonlinear  theory  to  describe  the  pinching  process  of  a  cylindrical  jet 
of  fluid  surrounded  by  a  second  fluid  with  possibly  diflFerent  density.  The  interface  between 
the  two  fluids  can  support  surface  tension.  The  theory  is  based  on  the  smallness  of  the 
jet  radius  to  a  typical  axial  wavelength  of  the  jet.  This  excludes  jets  which  try  to  roll  up 
and  such  behavior  would  appear  as  a  slope  singularity  in  our  model  evolution  equations. 
As  an  illustration  of  such  behavior,  consider  the  equations  (Bl),  (B2)  for  a  jet  in  air  in 
the  absence  of  surface  tension  (144  =  0).  Equation  (B2)  is  an  inviscid  Burgers  equation 
and  X  forms  a  shock  singularity  after  a  finite  time  (for  definiteness  we  are  considering  the 
periodic  problem  with  smooth  initial  conditions).  Similar  behavior  is  found  for  g-,  the  slender 
jet  approximation  is  violated  for  such  solutions  since  the  axial  length  scale  becomes  much 
smaller  than  the  radial  one  violating  the  original  assumptions  used  to  derive  the  model.  If 
surface  tension  is  included,  however,  the  equations  become  elliptic  and  allow  the  possibility 
of  admissible  pinching  solutions  described  in  this  article.  It  is  interesting  to  find  that  for  jets 
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surrounded  by  an  outer  fluid  so  that  the  Kelvin-Helmholtz  instability  mechanism  is  present, 
the  mode!  equations  even  in  the  absence  of  surface  tension  are  elliptic  and  pinching  solutions 
can  again  be  constructed. 

Finally,  we  return  to  the  similarity  solutions  constructed  in  Sections  5  and  6.  There 
appear  to  be  three  canonical  breakdown  regimes  corresponding  to  (i)  <  4/5,  (ii)  3  =  4/5, 

and  >  4/5  (see  Section  6).  In  what  follows  we  consider  the  validity  of  such  solutions  within 
the  context  of  the  long  wave  theory.  These  equations  are  derived  for  the  following  unsealed 
lengths  (stars  are  used  to  denote  unsealed  quantities); 

r*~I.  (7.1) 

As  a  finite-time  singularity  is  encountered,  it  is  seen  from  (4.1)  that  the  length  scales  over 
which  similarity  solutions  are  constructed  are: 

^  (if  T  r*~(i, -<)".  (7.2) 

We  know  from  (4.3)  that  a  =  4  —  4/1  which  when  substituted  into  (7.2)  along  with  the 
ordering  z*  »  r*  in  order  to  validate  the  slender  jet  approximation,  yields  the  ordering 

»  (is  -  (7.3) 

€ 

The  three  canonical  cases  (i)-(iii)  described  above  can  be  seen  in  terms  of  (7.3).  It  is  eas¬ 
ily  verified  that  cases  (i)  and  (ii)  satisfy  (7.3)  for  all  t  arbitrarily  close  to  tg.  Case  (iii), 
characterized  by  /9  >  4/5,  is  analyzed  by  letting  /i  =  4/5  -|-  6  with  <5  >  0;  this  requires  the 
ordering 

or  (7.4) 

Equation  (7.4),  therefore,  defines  a  range  of  times  for  which  convex  interfacial  shapes  are 
formed  at  break-down.  From  a  physical  point  of  view  we  expect  solutions  corresponding  to 
cases  (i)  and  (ii)  to  be  more  appropriate  when  the  ordering  (7.4)  is  not  observed.  Composite 
or  patched  solutions  are  possible,  however,  whereby  for  times  given  by  (7.4)  a  convex  drop 
begins  to  form  which  is  then  matched  to  a  concave  or  linear  shape  corresponding  to  fi  < 
4/5  and  =  4/5  respectively,  to  allow  evolution  well  into  the  break-down  regime  without 
violation  of  the  slender  jet  ansatz. 
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A  Derivation  of  the  energy  conserved  quantity 


The  flow  is  governed  by  the  incompressible  Euler  equations  which  are  conservative  and 
so  the  total  energy  of  the  system  is  expected  to  be  conserved  throughout  the  evolution. 
Contributions  to  the  energy  come  from  (i)  kinetic  energy,  (ii)  potential  energy  (in  the  presence 
of  gravity),  (iii)  interfacial  energy.  For  a  periodic  flow  of  period  d  (non-periodic  flows  simply 
alter  the  integration  limits)  the  total  energy  using  the  dimensional  variables  of  Section  2  is 

fR2  rlir  rd/l  J  fS  f'l-r  td/i  J 

s=  /  /  -p,\V4>-,\'rdrdedz+  /  /  -p,\Vcj>,\^drd0dz 

JS  Jo  J-d/'i  I  Jo  Jo  J-d/2  L 

-t-  /  /  /  gpirzdrdOdz  +11  I  gp\rzdrdddz 

Js  Jo  J-d/2  Jo  Jo  Jd/2 

r2jr  fdf2 

+(J  /  SdOdz.  (Al) 

Jo  J-d/2 

Clearly  the  potential  energy  contributions  are  zero  by  virtue  of  the  choice  of  reference  frame. 
Non-dimensionalizing  (Al)  according  to  (2.1)  (and  dropping  bars)  yields 


+  W,  S  dz.  {A2) 

J-i/2 

Making  use  of  expansions  (2.6a,b)  along  with  the  leading  order  solutions  (2.7)  and  the  change 
of  variables  (2.15)  yields  the  following  conserved  quantity 


£o=l^  (j/?(l  -  77)r/'V  + jT/(l  -  7/)'^X^  +  — 7^^^)  (A3) 

J-i/2  \4  4  a  J 

where  Zq  is  a  constant.  The  limits  of  integration  are  normalized  to  be  0  and  27r  as  they 
appear  in  equation  (3.2). 

We  have  shown,  but  for  brevity  do  not  include  here,  that  (A3)  is  consistent  with  the 
evolution  equations  (2.16),  (2.17)  for  the  equal  density  case.  This  is  achieved  by  taking 
the  time  derivative  of  (A3)  and  eliminating  r\i  and  from  (2.16),  (2.17)  to  show  that  the 
integrand  is  an  exact  ^-derivative  and  so  the  expression  is  identically  zero. 
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B  Equations  for  the  collapse  of  a  jet  in  air 


The  analysis  of  Section  2  is  used  with  a  zero  density  surrounding  fluid  which  implies  that 
(i)-2  =  0.  The  equations  to  be  solved,  to  leading  order,  are  (2.9)  and  (2.11)  with  B2  —  0.  If 
we  define 

^  =  *^o  X  =  Bx,, 

the  coupled  system  becomes 


Vt  +  {vx)z  =  0,  iB\) 

+  =  (B2) 

The  equations  are  in  conservative  form  and  so  the  integrals  (over  space)  of  7/  and  \  are 
conserved  quantities  (c.f.  (3.1a,b)).  The  energy  conserved  quantity  is  easily  obtained  from 
the  results  of  Appendix  A  and  in  particular  equation  (A2).  The  result  is 

o  ...  ,/A  .  _ 


^  Jo  ^  dz. 


We  note  that  (Bl)  and  (B2)  are  equivalent  to  equations  (2.1 1)  and  (2.10)  of  [TK]  lespectively 
if  (2.10)  is  dilTerentiated  with  respect  to  z  and  the  equations  are  put  in  terms  of  our  variables. 

As  a  check,  if  the  time  derivative  of  (B3)  is  taken  and  rjt,  Xt  are  eliminated  by  use  of 
(Bl),  (B2),  the  resulting  integral  is 

~\L 

which  is  zero,  as  expected,  since  it  is  an  exact  derivative. 

If  we  look  for  singular  solutions  of  the  form  (4.1),  the  orderings  (4.3)  emerge  as  well 
as  the  scaled  conserved  quantities  (4.4).  Mathematically,  therefore,  the  two  problems  are 
equivalent  and  the  same  restrictions  on  /3  apply.  The  equations  governing  collapse  in  the 
self-similar  region  are,  to  leading  order,  the  system  (4.6a,b). 
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Vortex  sheets  zero  surface  tension 

^=2/3  possible  g(f) 


FIGURE  1  Similarity  solutions.  Zero  surface  tension,  We  —  0.  Possible  solutions  for 
scaled  axial  velocity  jump  g(^)  for  ()  =  2/3. 


We*0  f  for  different  beta 


-5.0  -3.3  -1.7  0.0  1.7  3.3 

FIGURE  2  Similarity  solutions.  Zero  surface  tension,  We  =  0.  Scaled  interfacial  shape 
/(^)  for  different  ft  =  2/3,4/5,9/10. 
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We=0  g  for  different  beta 


FIGURE  3  Similarity  solutions.  Zero  surface  tension,  Wg  ~  0.  Scaled  axial  velocity 
jump  g{^)  for  different  (i  =  2/3,4/5,9/10. 


Wc=0  f,  g 


FIGURE  4  Similarity  solutions.  Zero  surface  tension,  =  0.  Detailed  view  of  the 
solutions  /(^)  and  g{^)  for  ft  -  7/8,  and  ^  >  0. 


We=0  log(f) ,  log(g) 


109(3)  -  stope»0.S6934  (exact  0.571429) 


log(»)  -  slope.-0.142142  (exact  -0.142857) 


3.9  4.1  4.3  4.5  4.7  4.9 

log(xi) 

FIGURE  5  Similarity  solutions.  Zero  surface  tension,  Wg  =  0.  Logarithmic  plots  to 
show  the  asymptotic  behavior  of  /(^)  and  ^i(^)  for  large  ^  and  /i  =  7/8.  Slopes  computed 
by  least  squares  included  on  the  figure. 

We/a=1  g  for  different  beta 


FIGURE  6  Similarity  solutions.  Non-zero  surface  tension,  We  =  1.  Scaled  axial  velocity 
jump  g{^)  for  different  /3  =  2/3,4/5,9/10. 


FIGURE  7  Similarity  solutions.  Non-zero  surface  tension,  W,  =  1.  Scaled  interfacial 
shape  /(4)  for  different  (a)  =  2/3,  (b)  ^  =  4/5,  (c)  ^  =  9/10. 
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log(xi) 


We/a*1  log  plots  beta=2/3 

xi  ->  +  infinity 


log(g) 


slope=  1.995 
(asymptotic  value  2.0) 


slope=-0.497 
(asymptotic  value  -0.5) 


(a) 


We/a=1  log  plots  beta=4/5 

xi  ->  +  infinity 


log{xi) 


8  Similarity  solutions.  Non-zero  surfacp  tension,  Wf.  =  1.  Logarithmic  plots  to 
confirm  asymptotic  states  for  large  Behavior  of  /(O  and  g{^)  for  different  li:  (a)  =  2/3. 

(b)  =  4/5,  (c)  fi  =  9/10.  Slopes  computed  by  least  squares  included  on  tlie  figure. 


24 


We/a=1  log  plots  beta=9/10 


log(xi) 


Figure  8c 
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